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^ ■ ABSTRACT 

^ , The strong lensing modelling of gravitational "rings" formed around massive galaxies 

is sensitive to the amplitude of the external shear and convergence produced by nearby 
mass condensations. In current wide field surveys, it is now possible to find out a large 
number of rings, typically 10 gravitational rings per square degree. We propose here, 
to systematically study gravitational rings around galaxy clusters to probe the cluster 
I mass profile beyond the cluster strong lensing regions. For cluster of galaxies with 

, multiple arc systems, we show that rings found at various distances from the cluster 

fSJ ' centre can improve the modelling by constraining the slope of the cluster mass profile. 

We outline the principle of the method with simple numerical simulations and we apply 
it to 3 rings discovered recently in Abell 1689. In particular, the lens modelling of the 
' 3 rings confirms that the cluster is bimodal, and favours a slope of the mass profile 

, steeper than isothermal at a cluster radius ~ 300 kpc. These results are compared with 

previous lens modelling of Abell 1689 including weak lensing analysis. Because of the 
difficulty arising from the complex mass distribution in Abell 1689, we argue that the 
, ring method will be better implemented on simpler and relaxed clusters. 
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1 INTRODUCTION 

In recent years, the modelling of cluster mass distribution 
has been progressively improved by (a) coupling the strong 
lensing (SL) analysis in cluster cores with weak lensing (WL) 
measurements at large radius (e.g Gavazzi 2005; Limousin 
et al. 2007b; Cacciato et al. 2006 ); (b) SZ measurements 
(e.g. Zaroubi et al. 2000 ; Dore et al. 2001; Seren 2007), (c) 
joint modelling of the cluster X-ray gas distribution (e.g. 
Mahdavi et al. 2007), (d) dynamical analysis of the veloc- 
ity distribution of the stars near the potential centre (e.g. 
Miralda-Escude 1995; Sand et al. 2004 ; Gavazzi 2005 and 
Koopmans et al. 2006). Despite these important improve- 
ments, it is not fully proved that cluster Dark Matter (DM) 
distribution closely follows the "universal" profile (Navarro, 
Frenk & White 1996, hereafter NEW) predicted by N-body 
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numerical simulations. Eor this profile, the DM space density 
is cuspy at the centre (r oc for r<rs) and behaves as 
outward. Measuring with great accuracy these two main 
characteristics is challenging and no consensus has arisen yet 
with the present day lensing observations. Recent analysis 
has shown that Sersic profile (Merritt et al. 2005) is also 
fitting the universal mass profile of numerical simulations. 
Importantly, the mass profile has only a weak dependence 
on the halo mass or cosmology, allowing to stack different 
measurements together to improve their significance. 

At large radius, analyses of the hot gas distribution from 
the X-ray observations (Pointecouteau 2005; Schmidt 2007), 
as well as weak lensing (Kneib et al. 2003) seem to favour 
an NEW-like profile for galaxy clusters. Eor elliptical galax- 
ies, Wilson et al. (2001) claimed that they were consistent 
with isothermal profile out to about 1 Mpc. Also, for mas- 
sive elliptical galaxies with gravitational rings, Gavazzi et al. 
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(2007) found that the WL slope of the mass density could 
be oc out to 300 kpc. 

In the core of mass concentrations, the actual existence 
of a DM cusp predicted by simple numerical simulations is 
more unclear. For spiral galaxy halos, the existence of a sin- 
gular density profile is still a debate because rotation curves 
are better explained by isothermal profiles with a core radius 
(Salucci 2003). Projection effects of non circular star orbits 
in triaxial halos has been invoked, but only in a few cases, 
to explain the linear increase of the velocity at the centres 
of galaxies (Hayashi et al. 2006 ). As a consequence, var- 
ious mechanical processes, such as gas cooling, supernova 
and AGN feedback, binary super massive black holes, dy- 
namical friction of the gas outflow during AGN activities, 
were investigated to explain the formation of a DM core ra- 
dius (Peirani et al. 2006 and references therein). For galaxy 
clusters, it is often argued that isothermal ellipsoid mass 
distributions with a flat core could better match the gravi- 
tational arcs geometry than NFW profiles (Sand et al. 2004; 
Gavazzi 2005; Gavazzi et al. 2003). X-ray observations gen- 
erally do not help much for this issue because, due to the 
limited spatial resolution of X-ray telescopes, one can just 
place upper limits on the radius of the smallest DM core, 
around 30 — 50 kpc (Chen et al. 2007). In summary, the DM 
density profile is still an open question, and high quality 
data are necessary to settle these questions. 

The deviation of light by masses is well described by 
gravitational lensing effect deduced from general relativity 
theory. The exquisite Hubble Space Telescope images (par- 
ticularly from the now defunct ACS camera) are providing 
the necessary lensing constraints to model observed gravita- 
tional lensing systems and may directly probe the existence 
of a "universal" density mass profile. However, one has to 
struggle to fully take into account the many observational 
parameters entering a lens modelling. 

Firstly, it is difficult to assess the stellar mass contribu- 
tion because the stellar mass-to-light ratio (M/L) is gener- 
ally badly determined, as well as the number of sub-halos 
and their galaxy occupation numbers (Wright et al. 2002). 
Secondly, lens modelling only probes the projected mass dis- 
tribution of lenses along the line of sight which introduces 
degeneracies in the 3D density profile if the mass profile is 
only determined on a small range of radius. When testing 
parametric models of mass distributions for a given DM con- 
densation, these current difficulties are only partly alleviated 
if we can probe the projected mass at many different radius. 
As an example, to disentangle between a flat core or a cusp 
with strong lensing it is not enough to detect and to analyse 
gravitational images very close to the centre, the so-called 
demagnified central images (Gavazzi et al. 2003) or inner ra- 
dial arcs (Mellier et al. 1993; Comerford et al. 2006). Even in 
such ideal cases, some information on the mass distribution 
beyond the Einstein radius is also critically needed. Beyond 
giant arc radii, one can use the information in the distortion 
of singly highly magnified arclets, in an intermediate shear 
regime (2 < /i < 3), also called flexion regime (Bacon et 
al. 2006; Massey et al. 2007). However, a generic difficulty 
similar to the one encountered in the weak lensing needs 
to be overcome. We do not know the shape of background 
sources and the flexion method must be used in a statisti- 
cal way. Only with the most deepest space-based observa- 
tions, it becomes possible to reach surface number density 



of background galaxies large enough to conduct such anal- 
ysis on a single cluster (see recent work on Abell 1689 by 
Leonard et al.(2007), where they reach a density of back- 
ground sources equal to ~ 200 sources/arcmin^). 

From the above discussion, we understand the difficulty 
to conduct an accurate measurement of the slope of the 3D 
density profile at large radial distances. Hence considering 
new probes of cluster mass profile is important. 

In this paper, we propose a method to investigate the 
slope of cluster mass profiles. Gravitational image systems, 
i.e. "rings" formed around galaxy cluster members, are used 
to analyse the slope of the cluster's density profile. Nowa- 
days, such rings can be systematically searched with ded- 
icated software (e.g. Gavazzi et al. 2007, in preparation; 
Cabanac et al. 2005). Here, we outline the method with 
three rings detected around Abell 1689. The main goal 
is to provide some constraints on the cluster potential at 
the location of the rings (i.e. at ~ 100" from the cen- 
tre of the brightest cluster galaxy). The coordinate sys- 
tem in this work is centred on the brightest cluster galaxy: 
a.720oo = 13 : 11 : 29.52, Jjaoo = -01 : 20 : 27.59. 

The paper is organised as follows. First we rapidly sum- 
marise the properties of gravitational rings observed in the 
field around elliptical lenses. In section [3] we illustrate the 
method by using simulated cluster profiles, lensing galaxies, 
and resulting images. These simulations show that we can 
put constraints on the local slope of the projected mass dis- 
tribution. In sectional the method is applied to Abell 1689. 
In this case we show that the three rings confirm that the 
cluster is dominated by a bimodal mass distribution and 
that the local slopes of both clumps are not much steeper 
than isothermal. Then, the results are discussed relatively 
to previous lensing models of Abell 1689, including a weak 
lensing analysis in the field of the rings. Finally, we conclude 
that this method should be better used on very relaxed clus- 
ters (single halo) with regular geometry to better probe the 
slope of the mass profile at various distances from the cluster 
centre. 

Throughout this paper we assume a cosmological model 
with r2m=0.3, f7A=0.7, Ho = 70km s"^ Mpc"^ At the red- 
shift of the cluster Abell 1689 (2 = 0.185) 1" is equivalent 
to 3.089 kpc. 



2 RINGS AROUND ELLIPTICAL GALAXIES 

The number of gravitational rings in the sky is large 
(Miralda-Escude & Lahar 1992; Blandford 2000). Observa- 
tions of strong lensing in the COSMOS field (Faure et 
al. 2007) have confirmed the estimation of their surface 
number density in optical survey (about ten per square 
degree with an average Einstein radius of about 0.8 — 
1.5"). With the development of large surveys like the 
SDSS (Sloan Digital Sky Survey) and more recently the 
CFHTLS (http:// www.cfht. hawan.ed u/Science /CFHLS/), 
several hundreds of such arc systems are currently being 
discovered (Bolton et al. 2006; Treu et al. 2005; Cabanac et 
al. 2005). Almost all of them are multiple images of back- 
ground galaxies found around massive ellipticals. Robot soft- 
wares are finding many of them among millions of objects 
in wide field optical surveys and can then be studied in de- 
tails with dedicated follow-up with HST or adaptive optic 
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systems (Marshall et al. 2007). Future space-borne surveys 
like JDEM/SNAP will allow us to discover as many as tens 
of thousands (Marshall 2005). 

Several structural properties have already been derived 
from the studies of elliptical lenses detected in the SDSS 
survey (Koopmans et al. 2006; Treu et al. 2006). It has 
been found that 1) the stellar mass is dominant within the 
Einstein radius and the velocity dispersion of stars almost 
matches the velocity dispersion of the lens model (a jai^ns = 
1.01±0.02 with 0.065 rms scatter); 2) the orientation of DM 
coincide with the light distribution within the Einstein ra- 
dius; and 3) the best fit for strong lens models coupled to 
a dynamical analysis of the star velocity dispersion demon- 
strates that an isothermal profile with p cx r-2±0 " can 
describe the total mass distribution at the Einstein radius. 

For elliptical galaxies in clusters, one can expect that 
the light distribution does not provide such a good geomet- 
rical description of the mass distribution. According to sim- 
ulation (Limousin et al., 2007c), the stellar contribution to 
the total mass is larger for galaxies in the cluster core. In- 
deed, as galaxies are crossing the higher density region of 
the cluster, the DM halo component is tidally stripped up 
to a radial distance which is not much larger than the galaxy 
optical size. This halo stripping has been studied observa- 
tionally by galaxy-galaxy lensing investigations. Although 
the deflection caused by galaxy scale mass concentrations is 
small (i.e. shear 7 ~ 0.01), it is measured in many clusters 
(Natarajan et al. 1998; Geiger & Schneider 1999; Natarajan 
et al. 2002a,b; Limousin et al. 2007a) and there is a clear 
evidence for truncation of galaxy dark matter halos in the 
higher density environments. The inferred average half mass 
radius is found around 40kpc in cluster cores, whereas half 
mass radii larger than 200 kpc are derived for field galaxies 
of equivalent luminosity (see Limousin et al. 2007a, for a re- 
view of galaxy-galaxy lensing studies). Hence, the DM mass 
and stellar mass can be of the same order within the ring. 
However in the following, we will only consider a total mass 
modelling of the ring. 

The recent SL2S survey (http:/ /www. cfht.hawaii.edu/ 
~cabanac/SL2S/) shows that the lens modelling of "ring" 
often requires the contribution of an external shear that are 
likely produced by nearby galaxy groups. In the outskirts of 
massive cluster it is also very likely that "rings" are present 
and are strongly affected by the cluster shear (Small et al. 
2007). For example in the HST/ACS images of AbeU 1689, 
several rings have been found (Limousin et al. 2007b). 



3 SIMULATIONS OF CLUSTER RINGS 

Cluster mass profiles can be described by different models 
such as the NFW mass profile, a truncated cored isothermal 
model (Kneib et al. 1996) or a cored power-law (Kneib et 
al. 2003; Broadhurst et al. 2005a). The latter two models 
are defined by 3 parameters whereas NFW is only defined 
by 2 parameters. For simple simulations we have chosen to 
describe the cluster mass profile by a single dark matter 
clump parameterised by a cored power law profile. The slope 
and core radius of the power law are chosen in such a way 
that it can produce giant arcs at ~ 45" similarly as in 
AbeU 1689. Ideally, the radius over which we will find rings 
should be significantly larger than the core radius so that 



the slope of the mass profile is almost constant at the radius 
of the "ring" . For a cored power law, the 3D density can be 
written as: 

p(r)=po(l + f )-" (1) 

so that the logarithmic slope of the 3D mass density p is 
given by n. 

Then, on top of this smooth cluster dark matter compo- 
nent we add a galaxy at a distance R from the cluster centre 
and we form a "ring" by lensing a distant galaxy. We then 
consider as constraints only the multiple images generated 
as part of the ring(s). In order to reproduce the astrometric 
accuracy of the HST-ACS images we add random errors for 
the position of the multiple images (±0.05" rms). In the fol- 
lowing, these artificial observations are then used to probe 
the slope n of the cluster potential without considering any 
additional giant arcs. Different observational configurations 
are investigated with one similar to the AbeU 1689 observa- 
tion discussed later. 

We will investigate the results on the measurement of 
n as a function of the ring distance to the cluster centre, 
and as a function of n. All the simulations were done using 
the LENSTOOL software (described in Julio et al. 2007) that 
allows to easily investigate using a Bayesian approach the 
important model parameters and their degeneracies. 

The rings in our simulations are described by a trun- 
cated isothermal profile for the total mass (PIEMD, Brain- 
erd, Blandford & Small 1996) 

'^^^ (l + rV^lo..)(l + rV^e\t) 

with an ellipticity of the total mass (e — (a} — l)^)/{a'^ + &^)) 
and position angle (PA) of the halo follow those of the light. 
In all simulations the cluster centre is set at [x = Q,y = 0) 
coordinates. 



3.1 Case A: Single Ring Configuration 

3.1.1 Optimising (n) 

To investigate with which accuracy we can recover the slope 
of the mass profile of the cluster, we have produced different 
mock configurations by varying the slope (tIin) and gener- 
ating multiple images around the lensing galaxy. Then we 
optimise the slope of the cluster mass profile using these mul- 
tiple images as constraints to estimate (nouT). With only one 
"ring", the maximum number of free parameters available 
is three for a triple image system and five for a quadruple 
image system. 

For the ring we can either choose a or rcut. In princi- 
ple, the velocity dispersion (Ji-ing can be measured from high 
resolution spectroscopy of the stellar component, so when 
there are only three parameters available, we thus choose to 
optimise rcut. Besides, the cluster slope n and the velocity 
dispersion Od of the PL potential 

4>{r,e) = 6^4^:^ro(l+(r/ro)2[H-ecos(2(e-eo))])('-">''' 

(3) 

are the parameters we will optimise, as the other ones (cen- 
tre, geometry, core radius) should be well known from the 
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n(in) 



Figure 1. Results for a single ring configuration. The ring is 
located 100" from the cluster centre. The deformations are gen- 
erated using a power law profile whose slope is nm reported on 
the horizontal axis. Using the constraints provided by the rings, 
we estimate this slope jiqut which is reported on the vertical axis. 
The black points with error bars correspond to a sample of sin- 
gle ring around a circular galaxy, red to a radial configuration, 
and blue to an orthoradial configuration. The black dashed line 
denotes the equality. 



modelling of the cluster multiple images. As our prelimi- 
nary goal is to see how well in this simplest simulation we 
can recover the input slope (tzh) of the cluster we present 
the analysis for three different geometric configurations of 
the lensing galaxy: (a) galaxy halo without any ellipticity; 
(b)galaxy halo with an ellipticity of 0.2 and PA=0° with 
respect to the line connecting the ring to the centre of the 
cluster (radial configuration) and (c) galaxy halo with an el- 
lipticity of 0.2 and PA varying up to 90° (orthoradial con- 
figuration). 

Results of the simulations and modelling are shown in 
Fig.[T] We have verified that the ellipticity and orientation of 
the lensing galaxy relative to the cluster shear do not change 
the accuracy in the recovering of the slope n. We plot the 
results for a fixed initial cluster velocity dispersion to show 
that the shallower the profile, the better the recovery of the 
slope. Indeed, the shallower the profile, the more massive the 
cluster, and thus the stronger its influence on a ring galaxy 
located at a fixed radius. 



3.1.2 Dependence with radial distance -Rring 

When optimising n, we have fixed the rings to be at 100" 
from the cluster centre because it corresponds to the radius 
of rings we have found in the field of Abell 1689. Finding 
rings further away from the critical region, where the shear 
of the cluster becomes weaker may bring interesting con- 
straints on the cluster mass profile. To investigate this pos- 
sibility with wider field survey of clusters, we have conducted 
simple simulation varying the ring radial distance from 100" 
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Figure 2. (Top): dependence of the accuracy of the method with 
respect to the distance between ring and cluster centre R(ring). 
In red, the input slope is n = 2.6, whereas in blue the input slope 
is n = 2.2. (Bottom): histogram of the probability density of n 
values for the point 250 " and n = 2.6 taken from the top figure. 



to 300" to see how well we can recover the local cluster mass 
profile. Results are shown in Fig. [J] The improvement in the 
determination of the slope does not depend on the azimuthal 
distribution of the ring around the cluster. But indeed if sev- 
eral rings are aligned on a same direction, we can not use 
them to estimate the cluster centre (see section [3.1. 3p . 

At distance larger than ~300" (~ 1 Mpc), as the cluster 
influence is getting small, it is only possible to give a lower 
limit on the slope n. Nonetheless, such a result can be use- 
ful as it would show that the slope of the mass profile has 
departed from the isothermal slope {n = 2). This can be un- 



derstood if we consider the total convergence k which is en- 
compassed within the ring radius R(ring). It conies from the 
contribution of the galaxy and the cluster (k — KG+Kdustor)- 
The convergence of the galaxy kg cannot be arbitrarily high, 
also the optimisation cannot accommodate any shallow clus- 
ter profile that will give a cluster contribution Kdustcr which 
is too high, since kg -I- Kdustor — 1- The better the likelihood 
on Kg is, the better will be the determination of n. Also 
the star velocity dispersion within the ring (total mass) can 
strongly improve the result. Only at large radius when both 
^duster and 7dustcr become very small variation across the 
ring, the method reaches its limit. (See Fig. [2} 

If a cluster has an NFW potential, the logarithmic slope 
of the NFW mass profile is "continuously" varying. A power 
law approximation with a core can match at the same time 
the total mass of the cluster within the ring and the local 
slope, but it results in an unphysical core radius rcore- Also 
for real clusters and dataset, the best way is to directly test 
an NFW model for the cluster component by including as 
many multiple images as one can find (see section 4). 

3.1.3 Finding the cluster centre 

In this section, we investigate how to retrieve the location of 
the cluster centre using several ring systems ~ 100" from the 
cluster centre. In this exercise, we take the cluster position 
as two free parameters. As shown in Fig. [31 a single ring 
cannot reliably find the centre of the cluster but does give 
some constraints on its direction. Note that the probability 
distribution contours are large and not centred on (0, 0). 

3.2 Case B: Multiple Rings Configuration 

We now consider three rings located at 100" from a cluster 
centre similarly as observed in Abell 1689. They surround 
the cluster: one is in the north of the cluster, one is south 
and the third one is west (see Fig. [S]). We then conduct the 
same simulations as done for a single ring but with the 3 
rings. 

Results on recovering the cluster potential slope are 
shown on Fig. O Compared to Fig. [T] we see that having 
three rings improves the measurement accuracy of the slope 
following a y/S factor. 

Regarding the cluster centre, we recover the centre of 
mass with the following accuracy x=0±50" and y=0±50". 
The constraints on the position of the clump are shown in 
Fig. [3] This measurement is not very accurate comparing to 
what we can achieve by analysing giant arcs systems gener- 
ated by a cluster. However, the simulation shows that not 
only are the rings sensitive to the slope of the potential, but 
also if we have enough rings, we can guess whether the loca- 
tion of the mass centre really dominated by a single massive 
clump. 

3.3 Conclusions from simulations 

A single ring can provide a good constraint on the slope of 
the cluster mass profile beyond its Einstein radius because 
the modelling constrains at the same time the local value of 
cluster convergence and shear. Schematically outlined, we 
can say that the cluster convergence enlarges the Einstein 
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Figure 3. Recovering the cluster position using the constraints 
provided by the ring(s). (Top): one ring (labelled by a blue star) 
is located at 100" from the centre of the cluster ( red star). (Bot- 
tom): similar results but with three rings (blue stars) located at 
~100". The contours define the Icr, 2cr and 3a region. 



radius of the ring and the shear rotates and changes the 
ellipticity of the critical line. Since both the convergence 
and the shear depend on the slope n at the location of the 
rings {R — Rring), the simultaneous adjustment of ^duster 
and 7ciuster provide a determination of n. The ring method 
can probe at which distance the slope of a cluster potential 
departs from n = 2 (isothermal slope). Using simulations, 
we are finding that there is a large range of cluster radii 
(-Re < R < lORe) over which rings can efficiently probe the 
mass profile slope. An optimum case would be to observe 
rings at various cluster-centric distances to better probe the 
overall DM profile. 

The ring method cannot put strong constraints on the 
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Figure 4. Results for a three ring configuration similar to 
Abell 1689. The rings are located 100" from the cluster centre and 
surround the main clump. The deformations are generated using 
a power law profile whose slope is rim reported on the horizontal 
axis. Using the constraints provided by the rings, we measure the 
slope nouT. 



position of the clump perturbation. The centre of the clump 
should be determined independently to take full advantage 
of the method (centre from X-ray or giant arc modelling). 
However, it can give information on the dumpiness of the 
DM distribution since each ring is more sensitive to its clos- 
est DM clump. With many rings surrounding a single domi- 
nant cluster potential, one can measure more accurately the 
slope of the mass profile, as well as determine the centre of 
mass. 



4 APPLICATION TO ABELL 1689 

With more than 31 arc systems, Abell 1689 (one of the rich- 
est clusters of galaxies at intermediate redshift) is produc- 
ing the largest number of strong lensing images. Although 
the X-ray map looks circular and is centred on the BCG, 
there is no cool core and the gas has complicated dynamics 
( Andersson et al. 2004) . Furthermore, the line of sight veloc- 
ity dispersion is complex (Girardi et al., 1997; Lokas et al. 
2006). Broadhurst et al. (2005b) found that the multiple arc 
systems can be modelled with a power law potential with 
n ~3.08, a surprisingly steep potential. As we will show be- 
low the 3 rings do not support a single power law potential. 

In fact, the most recent strong lensing models show a 
bimodal mass distribution (Miralda-Escude & Babul 1995; 
Halkola et al. 2006; Limousin et al. 2007b; Leonard et al, 
2007; et al. 2007). The cluster can be described as a domi- 
nant central dark matter mass clump (01 in the following) 
and a smaller perturbation associated with the north-east 
galaxy group (02 in the following). Although this two mass 
clump description is not fully satisfactory (see discussion in 
Limousin et al., 2007b), it can be considered as a first order 




Figure 5. Galaxy cluster Abell 1689. The field size is 280" X 280", 
corresponding to 1089 kpc X 1089 kpc. Gl, G2 and G3 are the three 
rings involved in this work. The images with lens galaxies sub- 
tracted are in the bottom row. The images we use as constraints 
are marked as little red circles. The diamond (C) is the optimised 
centre of the potential in the model we find with a single potential. 
It goes to the direction of the luminous clump in the north east 
(labelled by a cross 02). The centre found by X-ray is marked by 
a cross (Ol). 



approximation for the analysis presented here. Despite its 
complexity, we will apply the method presented above on 
Abell 1689. Furthermore, Abell 1689 has been extensively 
studied, allowing us to compare our results to former stud- 
ies as a check for the ring method. 

In the following analysis we will explore if the con- 
straints provided by these three rings can probe the bi- 
modality of the mass distribution, as well as the slope of 
each clump at the location of the rings (about 100" cor- 
responding to 300 kpc from the centre) without any strong 
a-priori assumptions coming from previous arc modelling. 
Next, we will use the results of previous modelling as input 
parameters to check if they remain compatible with the ring 
configuration. 
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4.1 Ring constraints 

Figure [5]shows the whole ACS field of view with the location 
of the three rings, as well as three panels corresponding to 
each ring. The three rings in Abell 1689 are at ~100" (~ 
300 kpc) from the cluster centre surrounding the two main 
luminous components (Table [T|. 

The three lensing galaxies are cluster members (based 
on their redshift) and have similar colour as the cluster red 
sequence. The lensed blue features around the ring G3 have 
a measured spectroscopic redshift of z=1.91 from LRIS at 
Keck (Richard et al., 2007). For the other lensed features 
around Gl and G2, we will assume for simplicity a redshift 
z—1. As a consequence, the velocity dispersion that we will 
derive has a loose meaning because of the assumption on the 
lensed redshift. The constraints we are considering here are 
only those derived from the blue lensed features detected 
around each ring (see Fig. [5)l . For a system with N multiple 
images, we get 2(A'^ — 1) constraints. Therefore, these give 
us a total of 16 constraints (Table [T|. 

In the following we only use these 16 constraints in or- 
der to probe the cluster without any a priori coming from 
the arc systems. The lenstool reduced quantifies how 
well we reproduce the location of the gravitational ring im- 
ages for each family of lens models. To describe each ring 
system, we use a PIEMD profile (Kassiola & Kovner, 1993; 
Limousin et al. 2005) for which we fix the position, ellipticity 
and position angle according to the light distribution. Then, 
we optimise the velocity dispersion a and cut-off radius Vcut 
for each ring. For ring G3, we also optimise the core radius 
Vcore to eusure a better match with this special lensing con- 
figuration although this has only a very minor impact on the 
results presented below. So in total, we use 7 free parame- 
ters to describe the three rings. This let us enough degrees 
of freedom to explore some cluster properties. 

4.2 Bimodality of the Cluster Mass Distribution 

If the central mass distribution of Abell 1689 is bimodal, 
the ring configuration will be different as compared to the 
one in a single central potential. The local external shear at 
each ring location has another direction. We have tried to 
simultaneously model the three rings with only one single 
central potential of either SIS, NFW or power- law profile 
for which we optimise the position of the centre and the 
dynamical parameters. We find that the P^r degree of 
freedom is always larger than 100! 

During this inconclusive optimisation, lenstool 
favours the centre of the potential to be between Ol and 
02 (Fig. [5]). This is unlikely since we know from both the 
lens model of the arcs and the X-ray map that most of the 
mass is centred on the BCG. 

Thus, we introduced a second mass component in the 
model. We fix all the parameters of the main mass compo- 
nent Ol to the values found in previous modelling (centred 
on the BCG) and we let free the location and the velocity 
dispersion of the second mass component. Much better solu- 
tions are found in this case, with a reduced smaller than 
4. In previous parametric strong lensing studies, the posi- 
tion of the second mass component is not well constrained. 
In this study we find a value (RA=13:11:31.360 ± 7("), y=- 
01:20:04.66 ± 21(")) close to the light centre of 02. This test 



confirms that the three rings are sensitive to the bimodality 
of the cluster central mass distribution, thus providing an 
extra evidence for bimodality in this cluster. 



4.3 The slope of the mass distribution 

4-3.1 Cluster modelling and rings 

The three rings are sensitive to the Abell 1689 gravitational 
field and in principle we can use them to constrain the lo- 
cal slope of the potential at the location of the rings. With 
the constraints of 3 rings, we cannot explore too many pa- 
rameters. Also we first assume that each mass component 

01 and 02 can be described by a power-law profile with a 
core radius. Thus, we fix the position of clump Ol to co- 
incide with its brightest galaxy, and the position of clump 

02 with the barycenter of its brightest galaxies. We also 
fix the ellipticity and position angle according to the clus- 
ter arcs modelling of Halkola et al. 2006 (table [2}. The only 
free parameters which are left free are the slope and veloc- 
ity dispersion for each mass clump (6 parameters for the 
two clumps). As discussed in section 4.1, we have 7 param- 
eters to describe the three rings. For the free parameters, 
the ranges of variation are chosen as follow. The parame- 
ters of each ring galaxy can vary between reasonable limits 
for a galaxy scale potential, i.e. a central velocity dispersion 
(Jo between lOOkms"^ and 300 km s"'^, and the cut-off ra- 
dius between 5 and 100 kpc. The velocity dispersion of each 
clump is allowed to vary up to 2500 km s~^. The slope n can 
vary between 1.0 and 3.0. Under these hypotheses the slope 
of the potential associated with the main central clump is 
found to be: n = 2.4lg4, which is a little steeper than an 
isothermal sphere. But surprisingly a large value of 108±32" 
is found for the core radius of 01. It corresponds almost ex- 
actly to the distance of the closest rings. Since we prefer to 
only use a power law approximation when the core radius is 
really smaller than the radial ring distance we try to force 
the core radius at the value of Halkola et al. 2006 (table [5}. 
Then we recover a similar value n = 2.46lE!;^g for the main 
clump but the mass of the second clump 02 becomes larger 
than the mass of 01. It appears that a larger core radius 
can match the mass found for each clumps with the mod- 
elling of the 31 gravitational arc systems but a small core 
radius for 02 does not. Clearly rings seem to tell us that 
the mass distribution of clump 02 cannot be approximated 
by a single central potential and that its mass distribution 
is probably more complex. This can be understood if clump 
02 is in fact a projection of a filamentary distribution of 
mass along the line of sight. We give as an illustration in 
Table [5] the results of the optimisation with the Halkola's 
core radius values fixed (in brackets). Besides, we will see 
below that we have a degeneracy of solutions between the 
velocity dispersion and r^ut for each galaxy. 



4-3.2 Degeneracies 

Parametric strong lensing models with many parameters ex- 
hibit various degeneracies. One reason is that strong lensing 
is sensitive to the projected enclosed mass which is degen- 
erated with respect to some of the parameters describing 
the 3D mass density. For a discussion on the degeneracies 
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Number RA DEC dBccC) Lens ^apec Lens magnitude i?e(kpc) Images 2 Number of constraints 

1 197.872 -1.309 114 0.1758 18.25±0.13 2.0±0.1 l.O(assumed) 6 

2 197.897 -1.345 90 0.1844 18.58±0.16 2.2±0.2 l.O(assumed) 4 

3 197.884 -1.369 110 0.1855 18.11±0.01 4.1±0.1 1.91 6 



Table 1. Strong galaxy-galaxy lensing events involved in this work: coordinates (J2000), distance to the BCG; spectroscopic redshift mea- 
surements when available; F775W AB magnitudes obtained from the surface brightness profile fitting by Halkola et al. 2006; circularised 
physical half light radius in units of kpc (Halkola et al. 2006); redshift of images; number of constraints provided by these images. 



encountered in parametric lens modelling, one can refer the 
work of Julio et al. (2007). 

In Fig. |6]we present the most relevant degeneracies for 
the slope of the main clump Ol and the other parameters 
found in this work. We see that the velocity dispersion and 
slope n of 01 are poorly constrained as expected. Similar 
degeneracies are obtained for 02. The bayesian framework 
in which lenstool proceeds, it is important to specify the 
meaning of the "best fit model", especially in such cases 
of strong departs from Gaussianity. Throughout this paper 
we define the best fit value for a given parameter as the 
mode of the marginalized distribution (integrated over all 
other dimensions). Hence we can see some unconsistancy 
between such values (summarised in Table and Fig. [S] 
which directly comes from bayesian techniques. For instance 
the high value for the slope 71(01) = 2.4lQ gQ is essentially 
due to the choice of the exploration range for (j(01). The 
inferred value for n{01) would be lower when decreasing the 
upper limit of the variation range of cr(Ol). For comparison, 
in a more frequent ist fashion, one would find values n{01) ~ 
1.5 near the minimum . Fig. |6]also shows the interest of 
measuring the velocity dispersion of Gl and clump 01 to 
break the degeneracy and get valuable constraints on the 
density slopes. 

We have noticed that the ring G2 is closest to 01. 
Therefore, we make a test only with G2 to see if it gives 
an equivalent result for the slope of 01. Like the single ring 
simulation in section 13.11 only the images around G2 are 
used in the optimisation. Now the slope and velocity disper- 
sion of the main clump (Ol) and the velocity dispersion of 
G2 are optimised with only 4 constraints. In this case, we 
take other parameters compatible with three ring modelling. 
A local slope is found to be equivalent to n — 2.4±0.2. For 
the second clump the degeneracies are larger and we found 
the logarithmic slope could be steeper than 2.4. 

Ultimately, one would like to try a simultaneous mod- 
elling of the 3 rings together with the 31 multiple arc sys- 
tems of Abell 1689. However, this would go beyond the 
scope of this work for two reasons. First, the description 
of the second clump is likely to be more complex than pre- 
viously thought. Second, a comprehensive investigation with 
LENSTOOL would need an enormous amount of computer 
time just to analyse one of several possible potential config- 
uration. Instead, we can test published models of Abell 1689 
to check how well the rings are reproduced (we call this the 
"ring test"). To do this, we took as an input the structural 
parameters of the cluster potential and we only solve for the 
ring configurations to evaluate the ■ As we have already 
demonstrated above, the power law model of Broadhurst et 
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Figure 6. Degeneracies between the slope of the main central 
clump and dispersion of clump one and Gl. The peak value found 
for the slope parameter is n = 2.46^g'gQ (see table These 
degeneracies underline the importance in obtaining the velocity 
dispersion of the clumps and lensed galaxies. 



al. (2005a) does not pass the test at all. For the Halkola et 
al. (2006) models, the ring test gives a reduced of 5.4 for 
the SIS model and 5.9 for the NFW model, suggesting that 
the local slope for both case shall not depart much from the 
isothermal solution. 
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main clumps Halo One Halo Two 

x(") [-2.1] [-33.1] 

y(") [-4.7] [56.7] 

e [0.159] [0.385] 

e [56.4] [106.3] 

Tcorei") [18.13] [24.28] 

cr (km/s) 709.15 ± 538.16 1119.97 ± 178.66 

rings Gl G2 G3 

x(") [3.527] [-86.088] [-38.581] 

y(") [113.703] [-14.64] [-102.812] 

e [0.492] [0.09] [0.443] 

e [-33.4] [-24.9] [43.8] 

rcorc (") [0.01] [0.01] 0.01 ± 0.07 

rcut (") 14.6±16.9 < 10 < 10 

cr (kms-i) 173.74± 13.99 162.07 ± 23.02 174.70 ± 27.40 

constraints 16 

free parameters 11 

X?ot 18.27 

X^'/d.o.f. 3.65 



Table 2. Best fitted parameters for the power-law mass distri- 
bution of the two clumps. The parameters in brackets are fixed 
during the optimisation. 



4.4 Conclusions for the Abell 1689 ring tests 

A preliminary conclusion is that modelling of the 3 rings 
seems to favour slope a bit steeper than the SIS value at a 
radial distance of 105", with n = 2.4lQ gQ. On the contrary, 
the best models of the multiple arc systems (current Einstein 
radius of about 50") favour an SIS or an NFW model with 
a scale radius rs below 40". Indeed, if is not significantly 
larger than the radius where most of the tangential arc(let)s 
are formed, the NFW slope remains close to n = 2 when r = 
rs thus very similar to an SIS model. Abell 1689 rings are at 
only two times the cluster Einstein radius and are thus in a 
situation where we do not expect a strong departure from an 
SIS slope. Hence the ring radius R(ring) is not large enough 
to derive stringent conclusion. Furthermore, the potential 
complexity of Abell 1689 is likely increasing the degeneracies 
of the solutions. In the following we will try to discuss these 
results and their compatibilities with a weak lensing analysis 
of Abell 1689. 



4.5 Discussions 

The lensed features around the three rings contain extra 
information on the mass distribution and mass profile of the 
cluster at a distance larger than the tangential arcs. Rings 
confirm the Abell 1689 bimodality, they also show that the 
slope of the mass distribution depart from isothermality at 
100" from the cluster centre. 

The three rings are located at i?ring ^ 105" from the 




Figure 7. The k map of Umetsu. The field size is 15' on a side. 
Blue contours are the reconstructed lensing convergence k. The 
lowest contour and the contour intervals are 0.05. The 3 rings 
used in this paper are labelled with red circles. Note that the 
local slope of k is remarkably compatible with the slope derived 
from the ring analysis. 



cluster centre. Considering that the dominant cluster mass 
profile can be characterised by the parameters of Ol only, 
which is legitimate given the difference in velocity dispersion 
between each component, we find the slope of the 3D mass 
density to be equal to — 2.4tp;gQ. Describing Abell 1689 with 
an NFW profile, this sets an upper limit on the scale radius 
rs of the NFW profile: rs < Rring- This upper limit on the 
scale radius can be compared with the weak lensing analysis 
of Abell 1689. Limousin et al., 2007b found both the strong 
and the weak lensing regime to agree in this cluster. From 
the CFHT data, they pursued a weak lensing analysis, and 
found a scale radius — 93". This value agrees with the 
upper limit provided by the ring constraints. We can even- 
tually try to go a bit further in this comparison between 
the ring results and the shear map of Abell 1689. Umetsu 
et al. (2007) have obtained deep images of Abell 1689 with 
the Subaru telescope and published a projected mass re- 
construction. From their figure 4 (See Fig [7]), we can see 
that their main centre is almost the centre C determined 
by the 3 rings(Fig [5]). Moreover, it is possible to estimate 
from their map the local slope n — 1 = — ALog(K)/ALog(i?) 
at the location of the 3 rings. We recover for Gl, G2 and 
G3 an equivalent value of n respectively equal to 1.95 ± 0.2, 
2.48±0.1 and 2.21±0.2. Although such a rough calculation 
of the slope is only valid for power law mass distribution, it 
seems compatible with the ring analysis. G2 is the closest 
ring to the centres of the two clumps and the slope obtained 
from G2 (2.4±0.2, see section l4.3.2|l is similar to the Umetsu 
et al. (2007) result. We can thus argue that a ring can bring 
similar information as weak-lensing analysis but at a fixed 
point. 

In summary, we found evidence for departure from 
isothermality at 100" from the cluster centre (Fig|6]), espe- 
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cially for the G2 location which is not perturbed by outside 
sub-clumps. Our results are not in disagreement with previ- 
ous lens modelling of Abell 1689. We found in particular that 
at the G2 location, the potential is steeper than isothermal. 

The three ring galaxies do have a similar F775W mag- 
nitude, between 18.11 and 18.58. It is interesting to notice 
that the modelling already infers comparable velocity dis- 
persions values, around 170 kms~^, and small cut-off radii. 
For rings G2 and G3, the cut-radii are found smaller than 
10". These central galaxies may be highly truncated, in good 
agreement with the tidal stripping scenario (Limousin et al. 
2007c, Limousin, Sommer-Larsen et al. 2007). To investi- 
gate further the exact value of cut-off radii, it is necessary 
to wait for good measurement of the stellar velocity disper- 
sion within the ring radius. 



5 CONCLUSIONS 

This work revisits the importance of external shear per- 
turbations for the modelling of galaxy lenses (see Dye et 
al. 2007). But instead of focusing our attention on the im- 
provement of the ring modelling by introducing an external 
shear, we probe the gravitational potential slope in the out- 
skirts of groups and clusters of galaxies thanks to the high 
sensitivity of the ring modelling both to the local external 
shear and convergence of the cluster. We test the princi- 
ple of the method with simple simulations using power law 
models to describe the cluster mass clumps. From this simu- 
lation we found that within the astrometric accuracy of ACS 
images a ring can be used to detect a logarithmic slope of 
the external mass density distribution up to about n = 2.8. 
For an NFW profile with a typical concentration parameter 
Cvir ~ 5 (Comerford & Natarajan 2007), this should allow to 
probe the cluster potential up to a large fraction of its virial 
radius. It is now possible to detect rings in the outskirts of 
clusters with new ground based surveys like CFHT-LS, but 
the determination of the direction and amplitude of the ex- 
ternal shear depends crucially on the high resolution HST 
images. On the other hand, the knowledge of the lens stellar 
velocity dispersion and of the redshift of the ring itself would 
increase the robustness to the ring modelling (Koopmans et 
al. 2006). 

Using the 3 rings found in Abell 1689 we have shown 
that it is possible to estimate the mass density slope. Hence, 
we have shown with the information given by rings alone 
that the potential of Abell 1689 is bimodal. With this lim- 
ited number of rings it is however not possible to better 
trace other DM substructures within each clump. Therefore 
a global modelling of the 31 arcs systems and rings remains 
a good challenge. 

We find a mass density slope a bit larger than n = 2 
for each DM clump at about 100" from their centre, but 
the complexity of the cluster potential really weakens the 
result. Since the most immediate appealing application of 
the method is to probe the slope departure from a mass 
density profile, we strongly suggest similar analysis on less 
complex clusters. The ideal lenses to probe the reality of a 
universal NFW profile would be a cluster or a fossil group 
with a bright cD galaxy having a single dominant halo, and 
which displays at least a distant ring, a multiple arc with 
known redshift, and a radial arc. Such a configuration would 



probe the potential at various distances from the centre. The 
discovery of such "golden lenses" requests a large survey of 
massive clusters. Most often clusters with a great number of 
arcs have multi-polar potentials (i.e. longer caustics) so that 
are more complex to analyse. However, rings are influenced 
by all their nearby environment (including foregrounds). In 
this sense, a statistical study of the environmental effect on 
the rings should be conducted in the future. 

Many rings can be discovered now in the field of deep 
wide field surveys. It will be possible to improve the method 
and to use it more systematically. As an example, the 
CFHT-LS "wide" survey will cover 170 degree , and we ex- 
pect to detect about 10 rings per square degree with an 
average redshift zicns = 0.65 (Cabanac et al. 2007). The 
clusters and groups in front of or at the same redshift as 
the lensing galaxies have approximately the same sky den- 
sity (Oguri 2006). Thus we can expect that several rings 
per square degree will be infiuenced by the external shear 
of a nearby mass condensation in the field. Such cases are 
already observed in the SL2S survey (Cabanac et al. 2007). 

For Abell 1689, we have found three rings within the 
small HST field (r ~ 100"), which may be consider as sur- 
prisingly high number. In fact, the large number of ellipticals 
in this cluster and the extra convergence that it adds to the 
lensing effect most probably compensate for the small field 
of view imaged with the ACS camer4i]- We also found two 
possible rings closer to the critical region of Abell 1689 {i.e. 
at radii smaller than 50"). However, we are not consider- 
ing these rings here, because these systems are in a crowded 
region and can not be described by a simple contribution 
of a galaxy and a cluster scale component. Any attempt of 
modelling these latter rings must involve all the multiple 
arc systems of Abell 1689. In wide field surveys, perturbing 
clusters may be found at a lower redshift than the ring- 
producing elliptical lenses (e.g. Small et al. 2007). In such 
cases, a proper analysis requires a multi-plan ring modelling. 

In conclusion, rings seem to be a promising tool to con- 
strain the mass distribution slope at large radii from the 
centres of groups and clusters provided dedicated observa- 
tions are carried out on a few golden lenses. They can pro- 
vide information on many structural parameters of halos of 
clusters and galaxies. In complement to the modelling of the 
multiple arcs in the cluster core, they can confirm or oth- 
erwise dispute the existence of the universal DM halo pro- 
file predicted by numerical simulation, as well as to study 
sub-halos and their cut-off radii. For these studies, it is cru- 
cial to measure the velocity dispersion of the lens galaxies 
which produce the gravitational ring images. It clearly ap- 
pears from the study of Abell 1689 that this method should 
be implemented first on a sample of very relaxed clusters 
or (fossil) groups in order to analyse a single DM potential 
with the simplest possible geometry. Only when large sample 
of rings is available, it will become possible to start a sys- 
tematic analysis of the external shear perturbation on ring 
shapes in correlation with more complex nearby (eventually 
foreground) environment. 



^ After the submission of this paper. King (2007) has shown with 
cluster simulation that the ring cross section is increased by a 
factor about 3 near critical lines. 
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